Towards a theory of attosecond transient recorder 
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Laser assisted photo-emission by a chirped subfemtosecond extreme ultraviolet (XUV) pulse is 
considered within an exactly solvable quantum-mechanical model. Special emphasis is given to the 
energy dependence of photoexcitation cross-section. The streaked spectra are analyzed within the 
classical picture of initial time-momentum distribution rj n j(p,t) of photoelectrons mapped to the 
final energy scale. The actual time-momentum distribution in the absence of the probe laser field is 
shown to be a poor choice for ri n i, and a more adequate ansatz is suggested. The semiclassical theory 
offers a simple practically useful approximation for streaked spectra. Its limitations for sufficiently 
long chirped XUV pulses are established. 

PACS numbers: 33.20.Xx, 33.60.+q 



O 
>Y 

^ : 



> 
00 

m 
in 

o 

OS 

o 



I. INTRODUCTION 

Subfemtosecond x-ray pulses are becoming a powerful 
tool for studying ultrafast processes in atoms [l|, [H, HJ 
and solids and their range of applications is growing 
rapidly Q . The characterization of ultrashort pulses is an 
important problem in this field. The attosecond metrol- 
ogy is based on the laser assisted XUV (extreme ultravi- 
olet) photoemission, whereby the temporal structure of 
the XUV pulse is reflected in the resulting streaked pho- 
toelectron spectrum [l|, Q • A subfemtosecond XUV pulse 
creates a photoelectron wave packet temporally confined 
within a fraction of the oscillation period Tl of the laser 
light. The wave packet is accelerated or decelerated by 
the superimposed laser field, and the dependence of the 
kinetic energy of the photoelectrons on the release time 
allows to record the time dependent spectra with a reso- 
lution around 100 as 

The gross features of the process can be understood 
within a quasi-classical model, in which the evolution 
of the photoelectron wave packet is split into two steps: 
first, a short XUV pulse creates an initial momentum 
distribution, which is then treated as a distribution of 
classical particles accelerated by the laser electric field. 
The model explains modulations of the width and cen- 
ter of gravity of the photoelectron spectra depending on 
the time delay of the laser pulse relative to the XUV 
pulse [l], 0. Quantum-mechanical calculations of laser 
assisted atomic photoionization performed in Refs. [1, 0] 
have confirmed that the XUV pulse duration can be 
measured by the width modulation. Recently, vari- 
ous quantum-mechanical formalisms have been proposed 
that reveal the connection between temporal character- 
istics of the pump pulse and the streaked spectrum and 
form a basis for characterizing both the pulse and the 
electronic system being excited [1, Q . Because the pump 
perturbation cannot be completely disentangled from the 
probe laser field the inverse problem - the reconstruction 
of the input pulse from the output spectra - is a difficult 



one. 

In Ref. [H an XUV pulse of duration 250 as was first 
measured by the streaked photoelectron spectra from 
neon atoms, and it has been shown that the streaked 
spectrum can be treated as a tomographic image of 
the initial time-momentum distribution of photoelectrons. 
The two-step model has become a common paradigm 
in laser assisted photoemission [l(| HH, U^|- Although 
the initial time-momentum distribution r- 1VL i(j>, t) appears 
to be an important concept in the theory of photoelec- 
tron streaking it is not straightforward to rigorously de- 
fine this notion. A quan tum-mechanical derivation of 
rini performed in Ref. [l3| showed that under certain as- 
sumptions a function independent of the streaking field 
can indeed be constructed, and it reduces to the time- 
momentum distribution ro(p,t) in the absence of the 
laser field. In Ref. flij an alternative prescription for 
r mi was introduced, which very satisfactorily described 
the streaked photoemission lineshape for a wide range of 
the XUV pulse durations and delay times. Understand- 
ing of the physical meaning and usage of the function 
nni(p, t) is important in order to correctly interpret the 
streaked spectra and obtain information about ultrafast 
electronic processes. 

In the present paper the classical approach based on 
nni(p, t), which hitherto has been used mainly for illustra- 
tive purposes, will be shown to have predictive power. We 
shall develop the semiclassical superposition approxima- 
tion (SSA), in which the pump and the probe pulses can 
be disentangled. For a bandwidth-limited pump pulse it 
reduces the calculation of the final momentum intensity 
distribution of photoelectrons J(p) to the integral [l4[ 



J(p) = 



+ 00 



U(t)I[ P -Ap(t)}S[ P -Ap(t)]dt, 



where the shape of the pump pulse enters through its 
time intensity envelope U (t) and spectral intensity enve- 
lope I(p), and the properties of the system enter through 
the excitation cross-section S(p). The effect of the laser 
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field is described by the momentum Ap(t) transferred to 
the electron created at the time moment t. The choice 
of the integrand will be explained in Sec. IIV A[ and in 
Sec. [Vjit will be extended to chirped pulses by replacing 
the function I(p) in the integrand by a function g(p, t) 
with an explicit time dependence: 



J{P) 



U(t)g[p-Ap(t),t]S[p-Ap(t)}dt. 



We shall consider a one-dimensional prototype of the 
photoexcitation of an atom or a localized state at a solid 
state surface and apply the SSA to link the input to the 
output using the initial time-momentum distribution as 
an intermediary. The paper is organized as follows. The 
computational method is described in Sec. [TIJ In Sec. IIIII 
we analyze the time-momentum distribution of photo- 
electrons from a chirped pulse in the absence of the laser 
field. Section [TV] discusses laser dressed photoemission 
by non-chirped pulses: it reveals problems in defining the 



initial time-momentum distribution function 



and 



troduces a practical ansatz for n n i. Section fVl generalizes 
SSA to chirped pulses and discusses the inverse problem. 



II. MODEL PARAMETERS AND METHOD 

In order to avoid any numerical approximations in solv- 
ing the time-dependent Schrodinger equation we employ 
the one-dimensional model adopted in Ref . ■ The un- 
perturbed Hamiltonian (in the units h = 1, m = 1/2) is 
H = —A + V(z). The model potential V(z) is a piece- 
wise constant function shown in Fig. [Ha) • The system 
is enclosed in a box, so that the spectrum of H is dis- 
crete. This setup describes the photoemission from a 
localized atomic state in a crystal close to the surface. 
We shall consider the lowest state of the model atom 
[ej = —18.1 eV, see Fig.QJb) and the figure caption]. The 
XUV pulse creates two wave packets traveling in oppo- 
site directions. We shall consider the spectrum of the 
one traveling to the right. Although the presence of the 
surface barrier between bulk and vacuum, see Fig. [JJa), 
is not essential for the present application, the barrier is 
retained for technical convenience: it reduces the total 
number of energy levels and slows down the right travel- 
ing packet. 

The electric field E(t) applied to the system is a super- 
position of the XUV (subscript X) and laser (L) pulses: 

E{t) =E x (t-t K )cQs[u x (t-t x )}+E h (t-t h )5m[u % (t-tJ]. 

The XUV and the laser pulse are confined to the intervals 
of length D x and D h centered at t x and t L , respectively, 
and their envelopes E x and E h have the same shape 



E x , h {t) = E : 
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see Fig.QJc). Note that D xx are twice the FWHM dura- 
tion often used in the literature. The intensity spectrum 
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FIG. 1: (a) Potential V(z) of the model ID system. The 
atom is represented by a 6 a.u. wide potential well centered 
at z = —8 a.u. relative to the surface barrier, (b) Density 
of states, (c) Superposition of the electric field of the XUV 
pulse of duration D x and the laser pulse of duration D L . 



I(oj) of such XUV pulse is close to Gaussian. To con- 
struct the classical model we shall need the normalized 
temporal intensity envelope of the XUV pulse U (t) de- 
fined by the formula 



U(t) = [E x {t)/El] 



(1) 



All calculations are performed for u x = 91 eV and 
E x = 10 6 V/cm, which is within the perturbational limit. 
We shall also consider linearly chirped pulses, with cu x 
changing in time: 



LO x (t) =LU X + b(t - t x ). 



(2) 



The parameters of the laser pulse applied to the system 
are typical of attosecond spectroscopy [l5j]: all calcula- 
tions are performed for frequency lo l — 1.65 eV (T L = 
2.5 fs) and duration D L = 5 fs. In the following we set the 
center of the laser pulse to zero, i L = 0. If not otherwise 
stated the laser field amplitude is E® — 6 x 10 7 V/cm. 

The perturbation in dipole approximation is zeE(t), 
and the time-dependent Schrodinger equation (TDSE) is 
solved in matrix form in terms of exact eigenfunctions of 
H. In spite of its simplicity the model is quite realistic. In 
particular, the energy variations of the excitation cross- 
section are in accord with those experimentally observed, 
see Ref. [l4j]. A more sophisticated development of the 
model has been recently applied to laser assisted photoe- 
mission from solids [16|. The wave packets created by 
the XUV pulse can be driven back to the ion and rescat- 
ter, thereby absorbing additional UV photons. This was 
studied for ionization in Ref. [TtJ and for electron scat- 
tering on an ion in Ref. [l8[ by solving the TDSE for a 
ID and 2D model system. In the present work we are in- 
terested in the energy distribution J(e) of photoelectrons 
recorded by a detector far away from the excitation re- 
gion. The spectrum J(e) is calculated by reexpanding the 
packet moving to the right in terms of the eigenfunctions 
of the system. 

We shall adopt the following convention about nota- 
tion. The three variables - energy e, momentum p, and 
frequency (XUV photon energy) u> - can be used one 
instead of another as argument of frequently used func- 
tions. They are simply connected: lu = e — ej, with 
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£i being the energy of the photoemission initial state; 
e = p 2 + Vq. The kinetic energy of the photoclcctron in 
vacuum is, thus, ei + uj — Vq. For the sake of brevity, 
we shall use e, p or u> depending on the context, e.g., 
J(w) is the short for J{u> + e{), and J(p) is the short for 
J(p 2 + V ). 

III. XUV EMISSION WITHOUT LASER FIELD 

The temporal evolution of the photoelectron spectrum 
during the XUV pulse is shown in Fig. [5] for a bandwidth 
limited and for linearly chirped pulses. The color map 
shows the time and energy dependent occupation number 
J(e,t) of eigenstates |e). The vertical cross section of 
the map at the end of the pulse J(e, t x + hD x ) = J(e) 
is the energy distribution curve (EDC) measured in the 
experiment. The time derivative of the function J(e,t) 
gives the energy resolved transition rate 

r (e,t)=dJ(e,t)/dt. (3) 

The function ro (p, t) is, thus, the time- momentum dis- 
tribution of photoelectrons in the absence of the laser 
field. 

In the absence of the laser field the observed final EDC 
obeys the simple formula [3]: 

J(e)~7(e)5(e), (4) 

where is the spectral intensity distribution of the 

XUV pulse, and S(e) is the photoemission cross-section, 

S(e) = \(e\z\e i )\ 2 . (5) 

Equation (fj| predicts the photoelectron spectrum with- 
out solving the TDSE. It is a consequence of the weak 
XUV perturbation and of the negligible probability of 
transitions within the (quasi-) continuum, e > 0, see Ap- 
pendix. (This relation is also obtained in the perturba- 
tional limit of the approximation to the TDSE introduced 
by Lewenstein et al. [19(, i.e., in the limit when the vector 
potential of the electric field can be neglected compared 
to the electron momentum.) For the present system it is 
found to hold for XUV pulses over a wide range of pulse 
parameters. Also the entire map J(e, t) calculated from 
Eq. (TJ| practically coincides with the exact result. 

The spectrum evolution is seen to strongly depend on 
the chirp rate. The effect of a linear chirp on the inten- 
sity spectrum can be characterized by the dimensionless 
index -q = bD x /T x , where T x is the full spectral width at 
half maximum (FWHM) of the non-chirped pulse with 
the same envelope and the carrier frequency lj x . When 
r\ approaches 4 the J(e, t) map starts acquiring a charac- 
teristic stripe structure, see Fig. [2fb,c,f). However, it is 
not straightforward to determine the value of the chirp 
rate from a known J(e, t) distribution. 

To illustrate the point, we show in Fig. [3] the time 
dependence of the center of gravity of the EDC e g (t) ob- 
tained by integrating along vertical lines in the transition 
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FIG. 2: (color online) Time dependent occupation of final 
states J(e, t) for the duration D x = 2000 as, V° = 3 eV (left 
column) and 400 as, V® = 15 eV (right column). The chirp 
rate b is zero in graphs (a) and (d), and 6 = 6 (b), 12 (c), 
96 (e), and 240 eV/fs (f). The photoelectron energy is given 
relative to ei + lu® [the frequency at t — is lu® = 91 eV, see 
Eq. {2J] and time relative to t%. The maximum intensity is 
set to 1. The vertical cross section at the right border of the 
map is the final intensity distribution J(e) recorded by the 
detector. 

rate map, see Fig. [5J The function e g (t) has a compli- 
cated form with variable slope, which depends not only 
on the photon energy sweep b but also on the pulse du- 
ration. The fact that b cannot be immediately inferred 
from de g (t)/dt is a combined effect of two factors: the 
FWHM of the EDC changes with time, and the broader 
the EDC the stronger it is affected by the energy de- 
pendence of the cross-section S(e). In the present exam- 
ple the cross-section is higher at lower energies, see right 
panel of Fig.O so the broader the frequency distribution 
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FIG. 3: (color online) Time dependence of the photoelec- 
tron spectrum center of gravity e g (t) for the chirp rate 6 = 
±24 eV/fs for three pulse durations: _D X = 400, 800, and 
1200 as. Positive (negative) chirp is shown by full (dashed) 
lines. Right panel shows the function 5(e), see Eq. ([5]). Pho- 
toelectron energy e is given relative to ei + lo® . 
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the larger the downward shift of the center of gravity. 
This produces a different effect on the spectra depending 
on the sign of the chirp, so the curves for opposite chirps 
in Fig. [3] are not symmetric. At t = D x /2 the curves of 
opposite chip converge, which follows from Eq. Q be- 
cause the spectrum I(oS) of a linearly chirped pulse does 
not depend on the sign of chirp. 

The function J(e,t) is the most fundamental quantity 
that contains information about the details of the ion- 
ization process. In the simple case considered above the 
properties of the system enter via the energy dependent 
cross-section S(e), and its structure is seen to be reflected 
in J(e,t). At the same time, this means that in gen- 
eral the XUV pulse cannot be satisfactorily characterized 
by gross features of the resulting photoelectron spectrum 
evolution. In addition, the function J(e, t) cannot be di- 
rectly measured in an experiment, which makes the prob- 
lem of disentangling the pulse characteristics from the 
properties of the system highly nontrivial. The only way 
to probe J(e, t) is the laser streaking technique, which we 
consider in the next section. 



IV. EFFECT OF LASER FIELD ON LINESHAPE 

Let us consider photoelectrons emitted parallel to the 
electric field of a linearly polarized light. They are accel- 
erated or decelerated depending on the sign of the mo- 
mentum transfer Ap(t) = (e/c)A(t) from the laser field to 
an electron created at the time t. As a result, the spec- 
trum shifts as a whole on the kinetic energy scale and 
distorts. One reason for the distortion is a finite width 
of the photoelectron spectrum: electrons with different 
momenta p acquire different energies, so the spectrum 
broadens when accelerated and narrows when deceler- 
ated. Another source of distortion is the finite duration 
D x of the XUV pulse: electrons created at different mo- 
ments t are differently accelerated according to the cur- 
rent value of the momentum transfer Ap(t). 



A. Initial time-momentum distribution 

In the absence of the laser field the momentum distri- 
bution of photoelectrons J(p) after the XUV pulse is over 
can be expressed as the integral of the transition rate: 



J(p) = 



r (p,t) dt. 



ho,. 



This formula can be generalized to the case of laser as- 
sisted photoemission by considering J(p) as a projection 
of the initial time-momentum distribution onto the scale 
of final momenta p, which is the essence of the classical 



picture of the streaking measurement Q: 

ix + ^-Dx 



J(p)= J r ini [p-Ap{t),t]dt. 



(6) 



This equation means that the electron created at a mo- 
ment t in the state p — Ap(t) is brought by the laser field 
to the final state p, and all such processes at different mo- 
ments t during the action of the XUV pulse are summed 
up incoherently. 

Let us concentrate on the physical meaning of the 
function r la i(p,t). This function characterizes the XUV 
pulse, and at the same time it depends on the proper- 
ties of the electronic system. To define this function 
requires a separation of the XUV and the laser pulses, 
which is, in general, not possible. For example, the time- 
momentum distribution in the absence of the laser field 
ro(p, t) cannot serve this purpose. First, even in the sim- 
plest case of a bandwidth limited pulse, for momenta far 
enough from the central momentum the function r (p, t) 
becomes negative because the momentum distribution of 
photoelectrons narrows with time. This is illustrated by 
the dot-dashed line in Fig. H^a), which shows two exam- 
ples of ro (p, t) calculated by the TDSE. The implication 
for the laser assisted photoemission is shown in Fig. 
which compares exact streaked spectra from the TDSE 
with those from the SSA. Although the contribution from 
negative values of ro is small it is not negligible, which 
is seen by the negative values of the spectral intensity 
in the dashed curves in Fig. [5] More important is that 
also around the central momentum, for which ro is ev- 
erywhere positive [dashed line in Fig. H[a)] , it provides a 
rather inaccurate final distribution J(p) (compare dotted 
and dashed lines in Fig. [5]). 

Instead of rp we can use in Eq. Q the separable ansatz 



r ini (p,t)^U(t)I(p)S(p) 



(7) 
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FIG. 4: (color online) (a) Temporal intensity profile of the 
XUV pulse envelope U(t), see Eq. (ft}, (solid line) compared 
to the laser- free transition rate functions ro(p,t) for the cen- 
tral momentum, po — £i + (jj x — Vo (dashed line), and for a mo- 
mentum far from the center pi (dot-dashed line), (b) Electric 
field _E L (t) of the laser pulse (dashed line) and the transferred 
momentum Ap(t). The circles indicate the release points t x 
of the XUV pulses studied in Fig. [S] 
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FIG. 5: (color online) Streaked photoelectron EDCs for XUV 
pulse durations D x = 1000 as (a, b) and 1600 as (c, d) and 
release points t x = 600 as (a, c) and -400 as (b, d) indicated in 
Fig.gfb). Dotted lines are the full solution of the TDSE, and 
solid lines are obtained from the SSA equation (fSj with the 
function n u i from Eq. (0. Dashed curves result from using 
the laser-free function ro in Eq. (]SJ. 



introduced in Ref. [l4| , which turns out to perform much 
better, see solid lines in Fig. \5\ In the phenomenological 
formula ([7J the explicit time dependence of the initial 
time-momentum distribution comes solely through the 
temporal intensity envelope U(t) of the XUV pulse [see 
Eq. (JTJ)] , which is shown by the solid line in Fig. 0Ja). 
(Oscillations with the frequency u x are, thus, completely 
ignored.) In the presence of the laser field the functions 
I and S in Eq. ^ also depend on time through the ar- 
gument p — Ap(t). The formula establishes a simple 
connection between the structure of the XUV pulse and 
the streaked photoelectron spectrum, and it can be useful 
for the characterization of the pulse provided it reason- 
ably simulates the true spectrum. It satisfies the two 
limiting cases: first, for an extremely short pulse, Ap(t) 
can be considered constant, and Eqs. © - yield a shift 
of the spectrum as a whole. Second, in the absence of the 
laser field the equation ^ is recovered. 

Figure Ufa) compares laser- free functions ro(e,t) with 
the function U(t) used in the ansatz ([7]). At the Fermi 
golden rule energy ei + lu x the true population profile ro 
is retarded by some 100 as with respect to the XUV in- 
tensity profile U(t), which is sufficient to cause an error 
of more than 3 eV in the location of the spectral max- 
imum. At the same time, equation ([6]) with the heuris- 
tic ansatz ([7]) reproduces the spectra with a surprisingly 
high quality (compare dotted and solid lines in Fig. [5]). 
Thus, in view of the sensitivity of the streaked spectra 
to small details of the function ri n i(e,i), the prescription 
(UJ appears to be a reasonable choice. The model per- 
forms well for durations well exceeding a linear interval 
of the Ap(t) function. It describes the dependence of the 
spectral width on release point and pulse duration, as 
well as the asymmetric shape of the streaked EDCs [see 
Fig. Eld)]. 



B. Effect of a rapid change of momentum transfer 

A detailed discussion of the laser field induced line- 
shape variations of photoemission spectra has been pre- 
sented in Ref. [HI ■ We shall now consider the possibility 
of inferring the XUV pulse duration from the measured 
streaked EDC. At a point where Ap(t) rapidly changes 
with time, electrons created at different moments t are 
differently accelerated according to the current value of 
Ap(t), so the longer the pulse is, the stronger the EDC 
broadens. This "temporal" broadening is accompanied 
by the energy stretch of the spectrum due to a uniform 
shift along the p axis (electrons with higher kinetic ener- 
gies acquire larger energy from the laser field than those 
with lower energies). If the integral momentum transfer 
can be neglected [around t — t L = ±840 as, see Fig. HJb)] 
the latter effect is minimal. It has been suggested in 
Ref. [H that in this case the temporal emission profile 
U(t) is uniquely mapped into a spectral distribution of 
photoelectrons J{e), and temporal information can be 
retrieved from a single streak record. 

Owing to the finite spectral width of the XUV pulse 
the streaking broadening Ar c (i.e., the difference be- 
tween the FWHM of the photoelectron spectra in the 
presence and in the absence of the laser field) depends 
on the pulse duration D x in a complicated manner. Fig- 
ures |UJa,b) show the dependence of Ar c on D x for two 
release points t x and three laser field intensities. Up to 
durations well exceeding the linear interval of Ap(t) the 
streaking broadening steadily increases with D x while the 
FWHM of the field-free spectrum steadily decreases. As 
a result, two different temporal profiles U(t) may give 
almost identical streaked spectra. 

The duration can, however, be immediately inferred 
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FIG. 6: (color online) Graphs (a) and (b) show the depen- 
dence of the streaking broadening AF C on the duration D x of 
the XUV pump pulse for two release points t x — 840 as (a) 
and 600 as (b) by the full solution of the TDSE (solid lines) 
and by the SSA (dashed lines). The curves are presented for 
three values of the laser field amplitude: Eg = 6 x 10 7 V/cm, 
as well as 2Eq and 4_Eo- Graph (c) shows the dependence of 
Ar c on release time for D x = 300, 400, 600, and 800 as. 
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from Ar e , and the SSA gives a reasonable approximation 
over a wide range of laser intensities. For small D x the 
width of the non-streaked spectrum increases, and the 
absolute values of Ap(t) start playing a more important 
role: the energy stretch (which varies in time) interferes 
with the temporal broadening. At the release point t x — 
600 as the integral momentum transfer is much larger, so 
the energy stretch at small D x leads to the presence of a 
minimum in the Ar o (D x ) curve. 

At short durations, the streaking broadening not only 
becomes small at the point of zero momentum transfer, 
but also the sensitivity of Ar e to the release point in- 
creases, as illustrated by Fig. HJc). This means that this 
point is less favorable for the measurement of D x than the 
measurement at the point of maximal momentum trans- 
fer. Figure [^c) suggests that Ar c measured at t x = 
(maximum Ap) and t x — ±1250 as (minimum), as well 
as the t x location of the extrema would characterize the 
XUV pulse more accurately. 

Figure M shows that the results by the SSA ©-© 
closely follow the exact solutions (except for the overes- 
timated maximal broadening for longer pulses), and their 
quality does not deteriorate with increasing the streaking 
field strength. 



V. EXTENSION OF SSA TO CHIRPED PULSES 

We now generalize the semiclassical approach to 
chirped pulses. We again rely on Eq. to separate 
the pump and probe pulses, but Eq. {?]) is not applicable 
now because it cannot discriminate between positive and 
negative chirp. Indeed, the pump pulse enters through its 
spectrum I(u>) and time envelope U(t), while both func- 
tions do not depend on the sign of the chirp rate b [defined 
in Eq. . In order to include the photon energy sweep 
in our model we must ascribe physical meaning to the 
formal notion of instantaneous frequency u> x {t). Thereby 
the spectral intensity distribution I{lo) in Eq. ([7]) will ac- 
quire an explicit time dependence. To achieve this we 
introduce an 'instantaneous' spectral intensity distribu- 
tion g(uj,t). In order to satisfy the basic equation in 
the absence of the laser field, this function should obey 
the relation 



g(oj,t)U(t)dt. 



(8) 



For a linear chirp and a temporal intensity distribution 
of Gaussian form, U(t) ~ exp[— (161n2) t 2 /D x ], the spec- 
tral intensity distribution I(lo) is again a Gaussian with 
FWHM T x = [64(ln2) 2 + D^b 2 ]^/D Xl and equation © 
has the explicit solution: 



g(uj, t) = exp- 



(oj - w° - btf 
7 2 /41n2 



2 h 2 
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_ r 2 Djb 

X A 



(9) 



Although in our case the function U(t) deviates from 
Gaussian we shall employ this form of g{u), t) and the 



requirement that the convolution ([5]) yield the FWHM 
r x of the actual spectral distribution I{oj). 

The semiclassical superposition (j6|) with the ansatz 



r iai (p,t) = U(t)g(p,t)S(p) 



(10) 



for the initial momentum distribution takes into account 
both the temporal and the spectral structure of the pump 
pulse. Thus, from the point of view of photoemission, 
the pump pulse is thought to produce at the time point 
t transitions with the spectral intensity distribution pro- 
portional to g(uj,t) [centered at u x (t)] and the maximal 
intensity proportional to U(t). In the limit of zero chirp, 
according to the requirement ([H]), the function g(u,t) be- 
comes time independent and reduces to I^lu), which leads 
back to Eq. ([7]). Note that in contrast to the pulselet 
representation of Ref. [8[, the function g(u>,t) does not 
contain phase, and its width is much smaller than the 
spectral width of a pulselet. 

Let us see how well the SSA © with the initial dis- 
tribution pTJj) describes the dependence of the streaking 
broadening Ar c on the chirp rate b of the pump pulse. 
For b = the broadening Ar o as a function of t x is prac- 
tically symmetric with respect to t x — 0, see Fig. Htc). 
For long pulses the "temporal" broadening mechanism 
dominates, and the curve Ar e (i x ) has two distinct max- 
ima (which merge into a broad one at t x = when the 
duration decreases). The symmetric shape is completely 

b = 10 eV/fs b = 25 eV/fs b = 50 eV/fs b = 100 eV/fs 




FIG. 7: (color online) Dependence of the streaking broadening 
Ar c of the photoelectron spectra on the release point t x for 
four values of chirp rate b = 10, 25, 50, and 100 eV/fs and 
four pulse durations D x = 400, 600, 800, and 1000 as. Solid 
lines are the exact solution of the TDSE and (red) dashed 
lines are obtained from the SSA equations ((6} and (|10|l . 
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distorted already at small chirp rates, see Fig. [7J One 
of the maxima is strongly damped and the other one is 
enhanced, so the shape of the curve provides a signature 
of the temporal structure of the XUV pulse. 

At a zero momentum transfer point the photoclcc- 
tron spectrum excited by a positively chirped pulse is 
stretched by the positive and compressed by the nega- 
tive laser field. This means that from the measured Ar c 
at the two t x points one can infer both the duration and 
the chirp of the pump pulse Indeed, Fig. [7] shows 
that both the broadening and the compression steadily 
increase with the pulse duration. However, the depen- 
dence on chirp rate is more complicated: for a fixed D x 
the amplitude of Ar o variations first increases with the 
chirp rate, and at larger b it starts decreasing. This non- 
monotonic behavior is caused by the interference of the 
energy stretch mechanism discussed in Sec. IIVBI over 
the whole XUV pulse duration electrons with a wide 
spread of energies are created, which are differently af- 
fected by the laser field depending on their energies. At 
higher chirp rates the spectral width T x of the pump 
pulse increases, which blurs the signature of its temporal 
structure. In addition, with increasing r x the energy de- 
pendence of the cross-section [factor S(p) in Eq. (fTTJ)) ] be- 
comes more and more important, whereby the laser-free 
photoelectron spectral width T° increasingly strongly de- 
viates from r x (see also a discussion in Ref. [3] and the 
comment on Fig. [3] in the end of Sec. IIIip . 

The SSA is found to well reproduce the spectral width 
variations for r\ = bD x /T x < 1.5, which covers a wide 
range of experimentally interesting pulses. For longer 
pulses or larger chirp rates it gives only a qualitative 
picture, see Fig. [7] The model systematically under- 
estimates the amplitude of the Ar c oscillations, with 
the largest discrepancy occurring when the energy sweep 
causes a compression of the spectrum. This has to be 
expected as the quantum interference between different 
processes that bring the excited electron to a given final 
state becomes most important in this case. The broad- 
ening of the spectra is reproduced much better; in par- 
ticular, the time point at which the spectral width is 
maximal is predicted rather accurately. The numerical 
experiments, thus, suggest that the limitations of the 
SSA lie in treating the interaction with the laser field, 
i.e., in Eq. ©, whereas the heuristic description of the 
energy sweep by the instantaneous intensity distribution 
in Eq. (flT)]) is quite reasonable. 

In view of the complicated interplay between differ- 
ent broadening mechanisms and because the pulse width 
r x is a non-linear function of both D x and b the question 
arises whether a chirped pulse can be characterized by ex- 
tracting only the gross parameter Ar e from the streaking 
measurements. A more general but complicated proce- 
dure is to measure a complete spectrogram J(e, t x ) and 
decipher it with a FROG algorithm (frequency-resolved 
optical gating), which under certain conditions allows to 
reconstruct both the intensity profile and the phase of the 
pulse [201 ] . A physically more transparent way is to estab- 
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FIG. 8: (color online) Dependence of the laser-free spectral 
FWHM r° and streaking broadening Ar e on the duration 
D x and chirp b. Positive Ar e correspond to the release point 
t x = —840 as and negative to t x = 840 as. The 6 = const lines 
(solid lines) are labelled by chirp rate in the units of eV/fs. 
The dots in the lines are placed at equal intervals of 50 as 
between D x = 300 and 1000 as. Points of equal duration are 
connected by dotted lines. 

lish a relation between gross features of the input pulse 
and the output spectrum and consider the release points 
at which the output is most sensitive to input. Follow- 
ing Ref. [3[, we shall consider zero momentum transfer 
points. 

In order to determine the temporal parameters b and 
D x from the spectral characteristics r° and Ar e the in- 
put values must be presented in the output coordinates. 
The map in Fig. [8] shows b = const lines parametrized 
by D x for the two zero momentum transfer points, see 
Fig. EJb). The data are obtained with the exact TDSE. 
According to our calculations, neither the b — const lines 
nor the D x = const ones intersect, which means that, 
in principle, both b and D x can be determined from a 
laser-free EDC and a single streaked measurement (for 
one of the Ap = points). The sensitivity of output to 
input strongly varies over the map: the uncertainty of D x 
increases with D Xl especially for small chirp rates, and 
at short durations the uncertainty of b increases. The 
complicated structure of the map implies that solution 
of the inverse problem requires a careful consideration 
of the photoexcitation process, as well as very accurate 
measurement of the EDCs. 



VI. SUMMARY 

The main result of this work is the extension of the 
well known factorization J(e) ~ I(e)S(e) of the photo- 
electron spectrum by a single short pulse to laser assisted 
photoemission [as expressed by Eqs. ([6]), |7J), and (jTUJ)] . 
This simple relation tells us that the laser free photo- 
electron spectrum is determined by the spectrum of the 
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exciting pulse but does not depend upon its temporal 
structure. In the streaking measurements it is just the 
temporal structure which is probed, and we have shown 
that the resulting spectrum depends on the intensity en- 
velope, while the phase structure of the pump pulse is 
blurred. In the case of a Fourier limited pulse it en- 
ters only indirectly, as the relation between the temporal 
and the spectral extent of the pulse. In the case of a 
chirped pulse it is insufficient to know the time envelope 
and the spectrum. We have shown that this more com- 
plicated temporal structure can to a good approximation 
be treated in spirit of the formal instantaneous frequency 
by introducing an instantaneous spectral distribution and 
again neglecting the actual phase relations. 

We have applied these considerations to calculating 
streaked spectra within the semiclassical superposition 
approximation. Although the concept of initial time- 
momentum distribution n n i (p, t) lacks a rigorous physical 
meaning a function with the required properties can be 
constructed, which is capable of describing the lineshapes 
of streaked spectra for chirped pump pulses. In addition 
to the temporal structure of the pulse it takes into ac- 
count its spectral extent, i.e., the fact that during the 
action of the pulse electrons are created in all states but 
with different rate. Both the temporal and the spectral 
factors are important for the shape of the streaked EDC, 
and the photon energy sweep, which both complicates the 
temporal structure and broadens the spectrum, provides 
a severe test of the SSA model. The model yields quan- 
titatively correct results for a wide range of parameters 
and provides a useful basis for the understanding of the 
role of pulse parameters and properties of the electronic 
system. The latter are included through the energy de- 
pendent photoemission cross-section, which was found to 
strongly affect the lineshapes. The property of the model 
to decouple the semiclassical calculation of the time evo- 
lution of the system from the quantum-mechanical cal- 
culation of the dipole matrix elements makes it attrac- 
tive for use for more complicated many-body systems, 
for which a consistent quantum-mechanical calculation 
is currently computationally infeasible: one can reason- 



ably suppose that the present formalism would remain 
valid with a more advanced description of the excitation, 
see, e.g. Ref. 2l| and references therein. 

The main limitations of the model naturally come from 
its neglect of the quantum interference that accompanies 
the acceleration of the wave packet by the laser field. In 
the case considered this led to underestimated broaden- 
ing and compression for sufficiently large pulse durations 
and chirp rates, but did not affect the qualitative picture. 

The phenomenological function r- ln i suggested in this 
paper differs from the time-momentum distribution r$ in 
the absence of the laser field in the two essential features: 
it discards the retardation of the final state population 
with respect to the pulse intensity profile and avoids neg- 
ative values of the population rate. Thereby the function 
rini provides the most direct connection between the in- 
put pulse and the output EDC, and, moreover, the two 
aspects are found instrumental in reproducing the shape 
of the streaked EDCs within the SSA. This means, how- 
ever, that in solving the inverse problem, i.e., obtaining 
information about the pump process from streaking mea- 
surements some essential features of the process may be 
lost. 

We have shown that in the simplest case of a linear 
chirp the relation between the input pulse and output 
spectrum can be presented on a map that establishes 
a one to one correspondence between gross features of 
the input and the output. This means that provided 
a reliable theory of photoemission is available and the 
experimental resolution is high the pump pulse can be 
characterized by a few streak records. 
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the diagonal ones. Let us retain in Eq. (|A.3[) only the ma- 
trix elements djj = (j\x\j) and djo — ( j | x \ } . Note 
that this approximation does not depend on the strength 
of the perturbing field E{t). Its validity depends only 
on the properties of the system. Then the equations for 
different j separate: 

ia 3 {t) = aQ(t)E{t)d jQ + aj{t)[tj + d n E(t)]. (A.4) 

Let us write aj{t) as a product cij(t) = exp[— iejt]£j(t). 
Then the equation for £j (t) is 



iij(t) = exp[iu3j t]^o(t)E(t)d j0 +^j(t)E{t)djj. 
Here ojjo = e 3 — eo- The solution of Eq. (|A.5[) is 
t 

Ci(t) = -i [ dt'Za(t')E(t')d j0 exp[w jQ t'] 



x exp 



idjj / E{t")dt h 



(A.5) 



(A.6) 



APPENDIX 



or, in terms of vector potential A(t), 



We shall show that the simple relation J(e) ~ I(e)S(e) 
is obtained in the perturbational limit eA(t) /c <C p of the 
time-dependent Schrodinger equation if one neglects the 
non-diagonal dipole matrix elements between final states. 
In order to avoid the necessity of treating the final states 
as plane waves we consider a bound system (e.g., a quasi- 
continuum as in Fig. [l} with the unperturbed Hamilto- 
nian Hq — ~A+V(x) and a time-dependent perturbation 
E(t) x. Let us write the Schrodinger equation 



iip(t)= H + E(t)x tp(t) 



(A.l) 



in terms of the discrete spectrum of the eigenstates | j ) of 
the unperturbed operator, Ho\j) = e j\j}- The solution 
consists in finding the expansion coefficients a,j (t) : 



£ j (t) = -i / dt' £ (t')E(t') d jQ exp[iuj ja t'} 



x exp[idjj{A(t') - A(t))]. 



(A.7) 



At this point we introduce the requirement that the per- 
turbation be small. Its implication is twofold: first, for 
A(t) — > the initial state can be considered unaffected by 
the electric field, so that £o(t') — 1. Second, the phase 
djj(A(t') — A(t)) can be neglected compared to Wjot'. 
Then reduces to the Fourier transform E(ui) of the 
XUV pulse multiplied by the dipole matrix element djo: 



£j(t) = -id j0 / dt'E(t') exp[iu j0 t']. 



(A. 



(A.2) 



The time-dependent Schrodinger equation then reads 

i$>i(*)l-?>= [Ho+E(t)x]j2^(t)\j)- (A.3) 

Let j = be the initial state, i.e., ao(0) = 1 and 
a,j(0) = for j / 0. It can be shown that for a suffi- 
ciently weak potential V(x) the non-diagonal matrix el- 
ements of the position operator x, are much smaller than 



The photoelectron spectrum after the pulse is over 
J(ujjo) = |£j(-D x /2)| 2 is then a product of the photoe- 
mission cross-section S((jJjq) — \djo\ 2 and the intensity 
spectrum /(wjo) = \E(uijo)\ 2 . (Of course, Eq. (|A.8|) can 
be applied to any t, i.e., it describes the evolution of 
the spectrum, too.) This result is also obtained in the 
perturbational limit of the equations of Ref. [l|| , but in 
the present derivation the assumption of plane-wave-like 
final states is recast as the approximation of vanishing 
non-diagonal elements of the perturbation between final 
states. 



